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Abstract 

The recent measurement of the third lepton mixing angle, ^13, has shown that, although small 
compared to 612 and 623, it is much larger than anticipated in schemes that generate Tri-Bi- 
Maximal (TBM) or Golden Ratio (GR) mixing. We develop a model-independent formalism 
for perturbations away from exact TBM or GR mixing in the neutrino sector. Each resulting 
perturbation scheme reflects an underlying symmetry structure and involves a single complex 
parameter. We show that such perturbations can readily fit the observed value of ^13, which is 
then correlated with a change in the other mixing angles. We also determine the implication for 
the lepton CP violating phases. For comparison we determine the predictions for Bi-Maximal 
mixing corrected by charged lepton mixing and we discuss the accuracy that will be needed to 
distinguish between the various schemes. 
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1 Introduction 



The observation of a non-zero value for the lepton mixing angle, ^13, challenges models of charged 
lepton and neutrino masses based on TBM or GR mixing. To date five reactor experiments, T2K[T], 
MINOS [2], DOUBLE CH00Z[3j, Daya Bay[4j and REN0[5j have reported measurements of the 
mixing angle. Taken together with the other neutrino oscillation experimental measurements two 
groups have recently produced global fits to the PMNS matrix mixing angles and Dirac CP violating 
phase, the most recent one is presented in Table [TJ 

Table 1: Results for the lepton mixing angles and Dirac CP violating phase taken from the global 
fit to neutrino oscillation data [6] (for previous fits see [8] ) . Separate results are shown for 623 for 
the cases of Normal Hierarchy (NH) and Inverted Hierarchy (IH) . 



Parameter 


Best fit Ifj range 3a range 


sin^ 6'i2/10-i 


3.0 2.87 - 3.13 2.7 -3.4 


sin^ 6'i3/10~2 


2.3 2.07 - 2.53 1.6 - 3.0 


sin^^as/lO"^ (NH) 
sin2 023/10-^ (IH) 


4.1 4.08 - 4.14 3.4 - 6.7 
5.9 5.68 - 6.11 3.4 - 6.6 


6/7r 


1.7 0.90 - 2.03 -2 



One may see that, at the 3a level, sin^ ^23 is in reasonable agreement with the TBM, GR and 
Bi-Maximal (BM) value of 1/2. The value of sin^ 612 is also in reasonable agreement at the 3a level 
with the TBM value of 1/3 or the two proposed GR values of 0.276 and 0.345 (but not with the BM 
value of 0.5) but the value of sin^ ^13 is not consistent with the TBM, GR or BM value of 0. 

However such models, based on a family symmetry, usually apply only to the contribution to the 
PMNS matrix from the neutrino mass matrix and there is an additional contribution from the mixing 
angle needed to diagonalise the charged lepton mass matrix. If the latter is given by an underlying 
SU{5) GUT, with the Georgi Jarlskog form [9j to accommodate the muon mass, one gets a contribu- 
tion to sin^ ^13 given sin^ ^13 = ^ ~ 0.5. 10~^, still inconsistent with the recent measurements at 
the 3a level. Thus, in order to fit ^13, the model must be modified to obtain a larger contribution from 
the charged lepton sector or the neutrino sector or both. Modifications of the charged lepton sector 
consistent with an underlying Grand Unified Theory (GUT) have been suggested |lllll21 [T3j by chang- 
ing the Clebsch Gordon factors associated with the GUT from those assumed by Georgi and Jarlskog. 
However a disadvantage of these schemes is that they lose the phenomenologically successful predic- 
tion for the Cabibbo angle mentioned above, the equally successful prediction |I4d/^s| = \/"W"7^ 
[14j as well as the prediction for the charged leptons Det = Det Mi that also follows from texture 
zeros. 

The alternative is to modify the neutrino sector. In this context one may note that TBM and GR 
models lead to a discrepancy in the small mixing angle, while providing a reasonable approximation 
to the large mixing angles. For this reason it is reasonable to consider small perturbations as the 

■^This is the GUT-related analogue of the phenomenologically successful prediction for the Cabibbo angle, sin 0c = 
\\/mdl'nia — e^^ yjm„/mc\ that results from a texture zero in the (1,1) entry and symmetry in the 12/21 entries of the 
up and down quark mass matrix |10) . 
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dominant origin of ^igQ. Our approach is to allow for mixing between the unperturbed neutrino mass 
eigenstates, concentrating on the bilinear mixing patterns that can keep small the deviation of the 
good TBM or GR predictions for 9i2. Such mixing often occurs in family symmetry models through 
higher dimension operators involving familon fields beyond the leading order. 

In this paper we will explore both the charged lepton and neutrino origins for ^13 in detail as well 
as the case where ^13 is partly driven by both sectors. This is done for the cases that the unperturbed 
structure corresponds to TBM, GR and bimaximal (BM) mixing. In Section [2] we briefly review two 
parameterisations of the PMNS matrix that will be useful in our analysis. Section [3] presents the 
TBM, GR and BM structures that provide the starting point for the perturbation analyses and 
Section m discusses the family symmetry structure of these models. Sections [5] and [6] discuss general 
neutrino mass perturbations about the symmetry limit that allow for non-zero ^13 for the case of 
approximate TBM, GR and BM mixing and derive the form of the Dirac CP violating phase and 
the mixing angles that result. In Section [7] we give our assumptions about the mixing angles that 
diagonalize the charged lepton mass matrix. Section [H] discusses how the results obtained still apply 
in cases that both the solar neutrino mass and the mass mixing (which are of similar magnitude) 
are generated as perturbations. In Section [9] we provide a quantitative estimate for the CP violating 
phase and constraints on mixing angles, derived using the fits to the remaining mixing angles at both 
the la and 3a levels of accuracy. This allows us to estimate the accuracy that will be needed to 
distinguish between the various cases. In Section [10] we discuss the case that the most general 3x3 
mass mixing preturbations are allowed and argue that, if the solar mixing angle is not to receive 
unacceptable large corrections without fine tuning, such cases are disfavoured and the simpler 2x2 
mixing case is preferred. Finally in Section [11] we present our Summary and Conclusions. 



The PMNS matrix arises from diagonalising the charged lepton and neutrino mass matrices, and 
involves the unitary matrices that transform the left-handed charged leptons and the left-handed 
neutrinos. These unitary matrices in general each contain 3 Euler angles and 6 phases and can be 
put in the form 



P is a diagonal phase matrix that can be immediately removed by rephasing the charged lepton 
fields. P'^ is a diagonal phase matrix (e*'''^ , e*'*'^ , e*'*'^ ) that can only be removed if the neutrinos are 
Dirac. Each matrix C/^j'^ is itself unitary and contains two parameters: a phase 6^'^ and a rotation 
by angle Ol''^ in the ij plane 



and acts trivially in the direction orthogonal to the ij plane. The resulting PMNS matrix is 



2 The PMNS matrix 



(1) 




(2) 



UPMNS = V'^V 



(3) 



^For recent reviews see [15M16| and extensive references therein. 
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which involves 6 Euler angles and 12 phases. 

On the other hand the PMNS matrix can be written in the PDG form 



UpMNS 



^ C12C13 S12C13 5136^'"^ 

-C23S12 - Sl3'S23Cl2e"^ C23C12 - Si3S23Sl2e"^ S23C13 \P, (4) 
\ S23'Sl2 - Sl3C23Cl2e'^ -S23C12 " Sl3C23Sl2e'^ C23C13 



where sia = sin0i3, C13 = cos0i3 with 6*13 being the reactor angle, S12 = sin0i2, C12 = cos0i2 with 
9i2 being the solar angle, S23 = sin 6*23, C23 = cos 6*23 with ^23 being the atmospheric angle, 6 is the 
(Dirac) CP violating phase which is in principle measurable in neutrino oscillation experiments, and 

ill -,§2 

P = diag(e ,e 2,0) contains additional (Majorana) CP violating phases /3i,/32. Note that this 
form involves only 6 physical parameters, {sij,5, (^2) and corresponding to the parameterisation 

UpMNS = R23U13R12P, (5) 

where Rij is a rotation in the ij plane by angle 9ij. 

In general the relation between the parameters of these two representations is cumbersome, but 
to leading order in sfj and Sj'3 it simplifies [171 IB] 

sl2^23e-^^'^^^+'-^ - ^'13^^36-''^^ 

«l2Cr2C^3e-^'- + si3Cr2S^3e^^'--'") 

43^^36-^'- (6) 



si3e" 




= Si3e 




si2e" 




= 8126 




S23e" 




= ■^236 





where 



'^13 = S + 

^23 = ^'2 

h2 = 0^-02 (7) 



and 



0x = /3i + (7i-73) 

^2 = /32 + (72 - 73). (8) 

Results for Sij are obtained by taking the absolute values of eqs ([6|), while 6 = S13 — 612 — 623- Note 
that the phases 7^ of P'^ lead to shifts in the Majorana phases /3i^2 but do not affect S or Sij. 



3 Models of lepton mixing 



In this paper we only consider models in which diagonalisation of the neutrino mass matrix gives 
bi-maximal atmospheric mixing and ^13 = at leading order. This class of models has a PMNS 
matrix of the form 



cos Q 


sin B 





— sin 


cos © 


1 


V2 


V2 


V2 


sin 


— cos © 


1 


V2 


V2 


V2 



-7, \P 



(9) 
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and the phase ambiguity discussed above allows us to adopt a particular phase convention in which 
all mixing angles are in the first quadrant. The associated Lagrangian for the neutrino masses is 
given by 

C = - {ma h'l + nn,h'l +mcvl) + h.c. (10) 

where the mass eigenstates are given by i'a,b,c = {na,b,c • which point in the particular directions 
ria b c in flavour space given by 



Ub = SQVe + Cq{v^ - Vr)l V2 

Vc = CQUe- Sq{v^- Ur)/V2. (11) 

3.1 Pure Tri-Bi-Maximal mixing 

Pure TBM mixing dSl [20] is generated by a PMNS matrix of the form of cqQ with tan 9 = 1 /\/2, 



I.e. 



Utbm 



1 1 1 



P = R23R12P (12) 



where Rij are the orthogonal matrices corresponding to the mixing angles, — l/^j ^23 = 1/2 
and S13 = 0. Assuming the charged lepton mass matrix is diagonal this form can arise from the 
neutrino mass matrix if there is a suitable non-Abelian family symmetry. In Section [5] we discuss 
how this matrix changes if one allows small perturbations in the neutrino sector about the TBM 
form. However in general we expect there also to be corrections to the TBM form coming from the 
diagonalisation of the charged lepton mass matrix and these must be included via cqs®. 

The structure of depends on the basis chosen for the states. In models of TBM mixing that 
follow from an underlying family symmetry one starts from the symmetry basis in which matter states 
belong to unmixed representations of the family group. In this basis we begin by taking = 1 so 
that the neutrino current eigenstates are i^e,ti,T: corrections from ^ 1 are discussed later. Then 
the pure TBM Lagrangian describing the neutrino masses has the form of eq ()10p where the mass 
eigenstates are given by eq (fTT]) with tan© = l/\/2, i.e. 

'2 





(f/x + 








{2Ve - 



+ z^r)/\/6. (13) 

In the absence of any perturbations we have pure TBM, with atmospheric, solar and the remaining 
orthogonal state given by i/@ = i/q, z/q = and = Uc respectively, with m@ = nia and m© = rrif,. 
The mass of the third state, m_|_ = nic, could be zero. The columns of Utbm are given by the 
direction vectors Uc^b^a- 

In the notation of Section [21 pure TBM corresponds to the case = 1 and C/fg = 1 in eq©. The 
combination f^23^i3^f2 '^^^ then be rewritten as P''i?23-Ri2-P'^', the phases in P'' can be absorbed in 
a redefinition of the charged lepton fields and those in P'^' can be absorbed in P'^. Thus one obtains 
V = R23R12P'' giving the TBM of eq(fT2]). 
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3.2 Golden Ratio mixing 

Another promising form for the CKM matrix that can be obtained from an underlying family sym- 
metry is the Golden Ratio case \21\ [22] given by eq([9]) with tan Q = l/cp where (j) is the Golden Ratio 
(1 + \/5)/2. This leads to the mixing angle di2 given by = '^/{V^'t') ~ 0.276, quite close to the 
fitted value. An alternative version of GR mixing[23l [23] has been proposed in which tan0 = 2/(j) 
corresponding to sin^ 9i2 = 0.345, again quite close to the fitted value. In both cases the mass 
eigenstates in eq pO]) are given by eq ifTT]) . 



3.3 Bi-maximal mixing 



For comparison we will also consider bi-maximal mixing |26t [271 ES] with the neutrino contribution 
to the CKM matrix given by eqQ with tan = 1 corresponding to 



BM 




(14) 



giving s\l 



S23 = 0.5 and S13 = 0, the reality of the matrix following from the underlying family 
symmetry. In this case an acceptable value for 612 requires large corrections from the charged lepton 
sector with, analogous to the quark sector, sizeable mixing in the (1,2) sector. Choosing 



BM 



one finds 



cos a 
e sma 




(15) 



Sl3 
^12 



'23 

5 



a/V2 

1/2 -|- a cos 

1/2 -aV4 
5K 



(16) 



4 Family symmetry structure 



In this section we discuss the discrete symmetry structure that leads to the unperturbed TBM, GR 
and BM schemes in the neutrino sector. For these cases the neutrino mass matrix has & Z2 ® Z2 
family symmetry. One Z2 is the same for all three schemes and is the -f-)- — r interchange symmetry, 
with generator 



U 



( 1 














(17) 



-1 



U has three eigenvectors given by the fa,b,c of eq (fTT|) with eigenvalues — 1,-|-1,-|-1 respectively 
and so the mass matrix invariant under U also has these three eigenvectors. The eigenvector with 
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negative eigenvalue is Va = (0, 1, l)/\/2 and leads to bimaximal atmospheric neutrino mixing in all 
three schemes. The other two eigenstates have positive eigenvalues and can mix with an undetermined 
angle 612, c.f. eqdH). Since Va is prevented from mixing with these two states, U is also responsible 
for setting 6*13 = 0. 

The remaining Z2 symmetry, with generator S, also has the states fa, and Vc of eq (|lip as 
eigenvectors: S Va,b,c = Pa,b,c^'a,b,c- The eigenvalues Ph and Pc have opposite signs, PbPc = — 1, 
preventing mixing between ui, and Vc- The general form for S for arbitrary is 



S = a /3 -7 (18) 




where 



a 



V2 



tan — cot 

{PgPb + 1) cot e + {PgPb - 1) tan e 
2(tane - cote) 



-{PgPb - 1) cot e - {PgPb + 1) tan 9 
^ 2(tane-cote) ' ^ ' 

A particular form for S, with a fixed value of 0, determines 012- 

Models built to generate these various mixing schemes are of two types (for reviews see I16j). 
The first type has a family symmetry, G/, that is broken to Gi, = Z2® Z2 in the neutrino sector. 
The second type has a family symmetry that does not contain S and U but has flavon fields, 4>g^b-, 
that, due to the family symmetry, acquire vevs pointing in the Ug^h directions, (l)g^h oc Ug^b- If these 
fields have opposite charge under a Z2 non-family-symmetry the mass eigenstates then are given by 
^a,b = 4'a,b ' ^- The Lagrangian then has the form given in eq (|10p and is invariant under S and 
U corresponding to an emergent Z2 ® Z2 symmetry. An advantage of this approach is that the 
underlying family symmetry readily applies to the quark sector too. 

4.1 Tri-Bi-Maximal mixing 

TBM can be generated by an underlying 5*4, which contains both S and [/, or by ^4 which contains S 
and, for a restricted non-family-symmetry breaking pattern, leads to a Lagrangian with an accidental 
symmetry generated by C^I^. For the TBM case S, for PgPi, = -|-1, has the form given by eq (jl8p with 
tane = I/V2 



p 

JTBM = 



/ -1 2 -2 \ 

2 2 1 
V -2 1 2 ) 



(20) 



Its eigenvectors i'a,b,c have eigenvalues Pg{+1,+1,—1) respectively. Thus, while the symmetry gen- 
erated by both U and Stbm is unbroken, there cannot be mixing between the states and the mass 
eigenstates are just Vgbc corresponding to TBM mixing. 



®See 1151 1161 and references therein. 
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4.2 Golden ratio mixing 



The first version of GR mixing, whicli we label GRl, has tanO = l/i?^ and can be obtained from an 
underlying family symmetry [30]; while the second version of GR mixing, which we label GR2, 
has tan Q = 2/4> and can be obtained from the dihedral group Diq [211 [25] . 

In the neutrino sector, the symmetry is broken to G^, = Z2® where the first is generated 
by the [i ^ t interchange generator \J of eq ()17p with eigenvalues +1,-1,-1. The second Z2 has the 
generator, S, given by eq IfTS]) . 

For GRl the generator with PaPh = 1 has a relatively simple form 



p 



/ -1 ^/2 -y/2 
\ -V2 1/(1) (t> 



leading to {Pa,Pb, Pc) = Pai+'i-, +1, -!)• The case with PaPb 

{Pa, Pb, Pc) =Pa{-l, +1,-1)- 



(21) 



-1 in eq (fT9l) is given by SgrU with 



4.3 Bi-maximal mixing 

Bi- maximal mixing can be generated by an underlying 54 symmetry [31] broken to = Z2 Z2- 
As above BM mixing in the atmospheric neutrino sector and the vanishing of 9^^ is driven by the 
first Z2 factor with the generator U of eq ()17p . The vanishing of the CP violating phase and maximal 
mixing for S12 is due to the second Z2 factor with generator given by eq ()18p with tan = 1 



1 



/ V2 -V2 

Sbm = 2 1 1 

\ -V2 1 1 

with eigenvectors i^afi^c given by eq ffTT]) with tan0 = 1 and eigenvalues +1,+1, 



(22) 



-1 respectively. 



5 Perturbations to TBM from the Neutrino Sector 

In TBM mixing the neutrino directions ya,b,c of eq (|13p play a special role and we therefore parame- 
terize deviations from pure TBM by allowing for small mixing between these states. There are three 
independent perturbations about eq (fT2|) arising from mass terms Uafb-, ^a^c and V},Vc and as a result 
the mass eigenstates are mixtures of the original eigenstates. In leading order this corresponds to 
the parameterisation given by 



-^(l + 3.* + 2t*) -i=(l-p-t) ^(1 + r-s) I (23) 
V ^(l-3s* + 2t*) --l.(l + |r*-t) 




where r, s and t are small complex constants that quantify the effect of VoVb^ VoVc and v^^Vc mixing 
respectively. 

These three perturbations correspond to a (small) breaking of the underlying Z2 Z2 symmetry 
to a Z2 subgroup generated by Stem, StbmU and U respectively. To generate a non-zero ^13 the 
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symmetry generated by U must be broken and so we do not consider the case of v^^c mixing. The 
residual symmetries Stem and StbmU prevent f^z^c mixing and thus keep small the deviation of 0i2 
from the tri- maximal form while allowing for a viable value of ^13. 

In an A4 theory I'a'^b mixing is natural in the sense that the U symmetry is an accidental symmetry 
and can be broken by an ^4 singlet family vev, leaving the symmetry Stem in ^4 unbroken. In 
5*4 the U symmetry can be spontaneously broken by an ^4 non-singlet family field but requires a 
mechanism to generate the appropriate vacuum alignment. In both A4 and /S'4 cases, the breaking 
that leaves a Z2 symmetry generated by StemU, corresponding to i^a^c mixing, requires a suitable 
vacuum alignment of a familon vev. 

For the case the symmetry is emergent the faZ/^ mixing is natural in the sense that it occurs if 
the Z2 non-family-symmetry, that prevents such mixing, is broken by higher dimension operators, 
perturbations to the vevs of (pa^b or from vevs of additional family-symmetry-singlet flavon fields. 
The generation of UhUc or I'a'^c mixing requires the addition of a 4>c familon vev. 

In the rest of this Section we assume that a single perturbation dominates. The cases Uai^b and Ua'^c 
mixing, equivalent to r 7^ and s ^ respectively, correspond to TMl and TM2 mixing discussed 
in [15] and lead to the so-called "atmospheric sum rules" relating phases and mixing angles that are 
discussed below. We discuss the case of more general mixing in section [10] where we demonstrate 
that, although it is possible to have more than a single mass-mixing operator, if the correction to 
tri-maximal mixing in the solar neutrino sector is to be naturally small, it is likely that a single 
mass-mixing operator should dominate 

5.1 UaUb mixing: Z2 Z2 Stem 

The Lagrangian describing neutrino mass is modified to have the form 

Lab = \ [M@{ua + VSenf + M^iub - V3e*uaf + M^vl) + h.c. (24) 

where e is a (small) complex expansion parameter. This case corresponds to the choice r = \/2e and 
s = f = in eq ()23p . Writing e = e"*"^' it is straightforward to construct Vy^. Comparing with 
eq([l]) one finds 



(523 ~ 2\/2si3sin(5' vr 

5^ + \/2si3 sin 5^ + tt 

^12 ~ — \/2si3 sinJ"^ 

^ „e 

Sl3 ''^ ^13 

§23 ^ \l/V2 + e sl^l 

s^a ~ I/V3. (25) 

Hence the Dirac neutrino phase is given by 6'^ = ^13 — ^23 ~ "^12 — Note that the initial unperturbed 
TBM form taken here is as given in eq (|12p and that this contains two unknown phases. This means 
that the Majorana phases are not determined and this conclusion also applies to all the mixing 
schemes considered in this paper. 
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5.2 i/qZ/c mixing: Z2® StbmU 
The modified Lagrangian now has the form 



(26) 



corresponding to the choice s = and r = t = in eq (j23p . Once again with the definition 

e = e~*^' sf3 and comparing with eq([T]) and extracting the Dirac phase one finds 



~ „<; 



^23 ~ |l/v2 — e * 8^3/21 
s^'s « 1/^3 

(5^ = (5^ (27) 
6 Perturbations to GR from the Neutrino Sector 

As for the TBM case a non-zero value for §13 requires that the Z2 factor generated by U be broken. 
To prevent large deviations of S12 from the GR value another Z2 factor should remain. For Va^b 
mixing PaPb = 1 in eq (fT9]l . which we call Sgr-, while for VaVc mixing PaPb = — 1 corresponding to 
SgrU. 

6.1 Uai^b mixing: Z2 -> Sgr 

The Lagrangian describing neutrino mass is given by 

Cab = \(M@{ya + ^^^b? + Mq{u^ - -^Uaf + M^vl) + h.C. (28) 

2 \ sm fc) sm H / 

Comparing with eq([T]) one finds 

^ „<: 
Sl3 ^ Si3 

s^3 « -L|i + e-'5^otesl3| 
Si2 ~ ■se 

= 6' (29) 
where, as always, we define e = e~''^^ s\^. 

6.2 z/qZ/c mixing: Z2 ^ Z2 ^ SgrU 
In this case we have 

Cac = ^ (Mmii^a + —i^cf + MqvI + M^{Vc + —Vaf) + h.C. (30) 

2 V ce ce / 
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giving 



6" 



Sl3 
1 



131 



(31) 



7 The charged lepton sector 

In order to compare TBM and GR perturbation predictions with experiment it is necessary to 
determine the contribution from the charged lepton sector. In this we are partly guided by the 
relations between down quark and charged lepton masses implied by an underlying GUT and we 
assume that the charged lepton mass matrix is hierarchical, with a similar structure to the quark 
mass matrices. 

Given the smallness of the equivalent quark mixing, we set = and drop the last term 
in eqs([B]). Similarly we choose S23 — O(^) analogy with the equivalent quark mixing angle, 
yC/CAf ^ 0{ms/mb). For the case of we will illustrate the possibilities by three choices. The 
first choice assumes the simplest possibility, namely = 0- The second guably the most 

plausible as it is the analogue of the successful relation in the down quark sector discussed above, 
is given by = \Jwt'' Finally, ^13 may arise entirely from the charged lepton contribution, with 

^ \/2 S13 |ll t ll2 1 [T3 l l29j. For example, this can result if the renormalisable (dimension 4) Yukawa 
couplings are forbidden by a symmetry and the coupling is generated by a higher dimension 5 term 
[33j. Although in this case the neutrino mass matrix is not perturbed away from TBM we include it 
for comparison with the perturbed cases. Note that although It is usual to compute 5 treating 
as a perturbation, in this case it introduces significant errors so in Section [S] we will keep the full 
expression. 



8 An alternative perturbative framework 

In the previous sections we have described the well-known TBM and GR mixing schemes, as well as 
their origin from family symmetries and perturbations about them that yield non-zero ^13. In these 
schemes there is no connection between the masses and mixing angles. Yet the observed value for 
^13 suggests an alternative scheme where only one neutrino mass arises at leading order. Let Va,b,c 
be some special directions in flavor space such that ^13 arises from ab (ac) mixing. For a normal 
hierarchy, the three largest neutrino mass operators are then 

= ^imal'a+'mbl'b) +^abJ^aT^b{'macJ^aT^c) + h.C. (32) 

The relative sizes of these mass terms are : nif, : niab {niac) = 1 : 0.16 : 0.27(0.19). Remarkably 
we see that the perturbation that induces mixing is larger than the solar mass term m^. This 
motivates searching for a new type of theory where only the atmospheric neutrino mass arrives at 
leading order, = maV^jl-, while the perturbation includes both solar and mixing terms, L'^ = 

rUb vll2 -h raab Va^b {rriac ^al^c)- 
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As a very simple example consider a theory where the neutrino fields v and two flavon fields ^a,h 
transform under a flavor group Gj, and there is an additional ^3 symmetry under which v and 
are singlets, but 0;, transforms as q, with = 1. In addition there is a flavor singlet field x with Z3 
transformation a that acquires a vev much less than the cutoff, giving a small dimensionless parameter 
e that transforms as a. The flavon fields acquire vevs of magnitudes Va^h in special directions n^,;, 
in flavor space, (^af) = Va,bn-a,b, that could be the directions of eqns (fT3|) for example. This leads to 
neutrinos i'a,b = {(pafi ■ z^) transforming as (l,a) under Z3. The Lagrangian for neutrino masses is 
then a perturbation series in the small field e = x/A 

Ci, = my Qci ul + e^C2 vl + e*C3 Va^b + ^le^)^ + ^-c (33) 

where rrii, is the overall neutrino mass scale, ci^2,3 are order unity dimensionless couplings, and A is 
the UV cutoff. The observed masses are accounted for by taking e ~ 0.2. In practice the predictions 
for the mixing angles and the phases are the same as for the cases discussed above but the scheme 
has the merit of relating the magnitude of the solar neutrino mass to the magnitude of ^13 in terms 
of a single expansion parameter e. This result would be spoiled if Va and Vb are not comparable, 
but could be regained by having nib and niab arise at the same order in Va and Vb, as occurs in the 
supersymmetric theory described by the superpotential 

y\^v = my[ -Ci I/„ + -C2 Ub + C3 VaVb + ••• I (34) 

where, for example, f and x are triplets under Gf and 4>a^b are anti-triplets, and the Z3 quantum 
numbers are as given above. 

In the next section we give our predictions for neutrino mixing angles and CP violation when 
ah and ac neutrino mixing is added to the TBM and GR schemes. A key point is that these pre- 
dictions are more general, resulting whenever neutrino masses are dominated by the three operators 
^a^a-ii^b^b-.i^a^bii^a^c)- They are independent of the organization of the perturbation theory, and 
depend only on the special directions Ua^b- 

9 Comparison with data 

We organise our analysis according to the size of s\2- 0, \/me/m^ or \/2si3. In all cases we include 
the error that results from a charged lepton mixing angle given by S23 = m^/m,-. 

9.1 Neutrino perturbations and = 

In this case the only correction from the charged lepton sector comes from 523- Prom eqs([6]) we have 

«13 = S?3 = «13 

(^13 = (^13 

5i2 = S12 

^23 = ^23 - 43sina/\/2 

•S23 = •523 + •523/2 + V2S23S23 COS a (35) 
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where a = b\'^ — is unknown. Thus 

6 = 5i3- 5i2 - 523 = + S23 sin a/ V2 (36) 
and the correlated uncertainties in 6'' and S23 are determined by 823- 

9.1.1 faZ^f, mixing 

TBM 

Prom eg ([25]) . 

6'^ = 6^ = cos-^( ^^'~°-^~''=^ ). (37) 
V2si3 

Using §23 from eq ()35p we determine the allowed range of the Dirac CP violating phase 5 given 
the allowed range of S23 from the fits of Table [1] and allowing for the uncertainty introduced by a. 
The results are shown in the first two lines of Table [2j 

GR 

Por the case of Golden ratio mixing eqs(|25p are replaced by eqs ()29p giving 

6- = 6^ = cos-i(^^tlll^^^l^). (38) 
2si3 cot B 

Using S23 from eq ([55|) to include the uncertainty introduced by S23, and taking tan 6 = (1, 2)/(/), one 
obtains the prediction for 6 in the (GRl, GR2) schemes shown in Table [2J 

9.1.2 fafc mixing 
TBM 

The difference compared to the case of faffc mixing is that eqs (f25]l are replaced by eqs ([27|) . In 
practice the only change in the analysis is that eq ([37|l is changed to 



^/2(4|-0.5-sy4), 



= = cos-^(-^-^^^^^^^ "^-i'-' ). (39) 

Sl3 



The resulting prediction for S is given in Table [2l 
GR 

In this case the only change is that from eq (|3ip we have 

6'^ = 6^ = cos-i(-^^t^^ZLAt^)_ (40) 
^ 2si3tane ^ ^ ' 

9.2 Neutrino perturbations and s'^2 = 

A non-zero affects the determination of s^^ and si2- Prom eq® we have 



u2 
'13 



Sl3 + ^ll/^ + V2si3s[2 COS f3 



613 = Si3 + s[2sm/3/{V2. 



Sl3) 



(41) 
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where (3 = + ^12 ~ "^13 • Using these equations for given S13 and may determine and the 

uncertainty in the relation between and (and hence the error in the relation 5 = 5^) in terms 
of the unknown phase j3 and thus estimate their (correlated) errors. 
Turning to S12, eq® gives 

VSsiae-*'^!^ = e-'^'i^ - A^e''^'^^ . (42) 

We see that the prediction for S12 depends on unknown phases so now there is a range of allowed 
values as shown in Table [TJ Also, due to the unknown lepton phase, there is now an uncertainty in 
the relation between the neutrino contribution to the phase and the full phase, 812 = (5i2 i s^^2) 
this translates to an uncertainty in the determination of 5 of about 2^ . There is also an error coming 
from the unknown lepton phase, a, as discussed in Section [9.11 

Having allowed for these charged lepton contributions the analysis proceeds as detailed in Section 
19.11 The predictions for the phases and angles are given in Table [2j 



9.3 No neutrino perturbations and — V^-^is 

In this case ^13 is entirely given by the charged lepton sector. Given the size of it is necessary to 
keep the full dependence on it rather than use the approximate eqs([6]) giving (up to the correction 
from the uncertainty in S23) 



|S23S12 - Sl3C23Cl2e"^| = 

S23C13 = (43) 

The result for S23 is shown in Table [5] where the error is dominated by the uncertainty in §23 
which we estimate as m^/m,-. 

9.3.1 TBM 

Currently the largest uncertainty in the above prediction for 5 arises from the experimental uncer- 
tainties in S12 and S13, giving a range ±(0.54 — 0.61)7r at la and ±(0.48 — 0.68)7r at 3<t, as shown in 
Table [2j Thus a pure TBM neutrino mass matrix together with a hierarchical charged lepton mass 
matrix leads to near maximal Dirac CP violation. The experimental uncertainty in 8 is currently 
dominated by the experimental uncertainty in 5^2- la the fractional error in 5^2 is about ±0.06 
and this leads to an uncertainty of ±0.037r in 8. A reduction in this uncertainty by a factor 2 would 
lead to a prediction of 8 at the 5% level. 



9.3.2 GR mixing 

We find the range ±(0.41 - 0.47)7r at la and ±(0.32 - 0.51)7r at ^a for GRl and ±(0.56 - 0.64)7r at 
la and ±(0.51 - 0.72)7r at J.a for GR2. 
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9.4 BM mixing 

From eq ()16p we may determine the mixing angles and Dirac phase. For S23 the error is dominated by 
the leptonic contribution and, for S23 = fn^/mr, this gives the range = 0.44 — 0.56. BM mixing 
is inconsistent with the allowed range of at the la level. At the 3a" level we find 5 = (0.87 — 1)tt. 



10 More general mixing 

Up to now we have only considered perturbations about pure TBM where a single mass-mixing 
operator, Uafb^i^b^c or fa^c, dominates. Of course more general mixing is possible and we illustrate 
this with two examples. 

10.1 A model based on the flavor group ^4 

Pure TBM results in a model based on the flavor group ^4 having three flavon fields (p, (ps and (px 
and leading interactions [34J 

C = j (^VeiMe^hd + y^TlYl^^ha + yr{<pTl)"T^hd^ + ^ (^Xa(t>{ll) + Xb{(t>sll)) KK + h.c. (44) 

The dominant perturbations to pure TBM result purely from the neutrino sector, and arise from 
three operators [M] 

^ {xc{<PT<Psy{ll)" + Xd{(t>T<Ps)"m' + Xe(t>{<pTll)) KK + h.C. (45) 

We find that the couplings Xc and x^ lead to only the Vai^c mass operator, while the Xe couplings 
leads only to Vb^c mass mixing. Hence in this model ^13 must arise from Xc or x^- If the deviation 
of S12 from tri-maximal mixing is to be small as required by the data, the coefficient Xe must be 
small. In this case the model reduces to a single mass-mixing operator. Dominance by more than 
one mass-mixing operator requires a correlation between the coefficients of the mass operators, as 
the next example demonstrates. 

10.2 A model that preserves = 1/2 and = 1/3 

Maintaining the tri-maximal solar angle and the maximal atmospheric angle requires 

\Ue2? _ 1 |C/^3P _ , 



|r/2 12 2' \JP |2 

I'^ell ^ I'^tSI 



(46) 



and occurs in the scheme |35 



Cm = {va + -^(z^b + \/2z^c)^ + ^Mq [ub - + " \/f ^ ^« j + ^-c- (^7) 

with simultaneous mixing of Vb and Vc with Va- This evades the correlation discussed above between 
the departure of S23 from TBM mixing and the value of siglf]. Although this requires a strong 

drawback of such a scheme is that such a modification of the familon vevs changes the leading order structure 
of the Dirac mass matrices of the charged leptons and quarks and thus spoils the phenomenologically successful mass 
predictions that follow from a (1,1) texture zero. The bilinear mixing discussed in Section [S] arises from higher order 
corrections to the Majorana mass matrix and thus does not spoil the (1,1) texture zero relations [36j . 
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correlation between the mixing of these states it has been argued that this can happen quite naturahy 
in an A4 model through vacuum alignment of the familon vevs [35]. To see how this may come about 
note that the structure of eq (j52|) follows from a modification of the structure discussed in Section 
14.11 namely 

Cm = ^Mc^ii^iOif + ^MQ{ui0lf + h.c. (48) 

where now the triplets vevs have the form (pa oc (e, 1,— 1) and (pb oc (1,1,1). The modification 
compared to the tri-bi-maximal mixing scheme is the appearance of the entry proportional to e in 
the vev of the familon field (pa- 

In the proposed alignment scheme there are additional triplet familon fields (pi, (ps and 1^23- AH 
the familon fields have vevs driven by radiative breaking through a potential of the form V = m?\(p\'^ 
where rm? becomes negative at some high scale through radiative corrections. This term is SU (3) 
symmetric so to determine the vacuum alignment one must look for potential terms splitting the 
degeneracy. An ^4 invariant that does this has the form |(/>"f*(/)jp. If its coefficient is positive 
(negative) the preferred vev is oc (0,0,1) (oc (1,1,1)). It was argued that this mechanism readily 
leads to the vevs 

I <</.i23| >oc(l,l,l), < |0i| >a (1,0,0), < |03| >oc (0,0,1). (49) 
Finally the alignment terms 

Al23|</'I23<^23P + Xl\<p\hz? + Al|</>I</'a|^ + \2-i\4>2Z<l'a? (50) 

were added to try to obtain the desired alignment. The first two terms with positive coefficients 
force < \(p2'i\ >ot (0, 1,1). The last term forces (p22, and (pa to be orthogonal but does not require the 
second and third terms of ^23 to be non-zero. This is determined by the third term - if Ai is positive 
the first term vanishes while if it is negative the second and third terms vanish. Thus these terms do 
not drive the desired form of the vacuum alignment, (pa oc (e, 1, — 1). To arrange for such alignment 
requires a modification of the alignment terms. The simplest possibility arises if (pi and (p^ have the 
same quantum numbers under the symmetries beyond that are usually introduced to limit the 
allowed form of the Lagrangian. In this case the second and third terms of eq (|50|) can take the forrG0 

All (01 + X(P^)^2Z? + All (01 + y^z)Ua? (51) 

where x and y are constants. Now it is straightforward to see that, for small y, the vev of (pa has 
the desired form with e = —y < (pi > / < (p^ >. However to avoid spoiling the alignment of (P23 it is 
necessary that x be large and to keep the correcting within the oserved limits requires y/x < 10~^. 
Given that one expects the coefficients x and y to be of 0(1) this looks an unnatural requirement. 

In summary, while it is possible to achieve the vacuum alignment necessary to achieve the form of 
eq ()50p . this example shows that it requires a very complicated alignment mechanism and even with 
this requires some fine tuning of 0(1) coefficients. Thus, although it is possible to have more than a 
single mass-mixing operator, if the correction to tri-maximal mixing in the solar neutrino sector is to 
be naturally small, it is likely that a single mass-mixing operator should dominate and the structure 
analysed in the previous Sections apply. 

^The most general structure is more complicated but these terms suffice to illustrate the point. 
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11 Summary and conclusions 



In this paper we have explored the possibihty that the structure of the neutrino mixing matrix is 
given in zeroth order by the most promising family symmetries that have been suggested, namely 
those leading to TBM, GR or BM mixing and perturbed by mixing between the zeroth order mass 
eigenstates, fa, i^b and Uc- 

We discussed the underlying Z2 x Z2 family symmetries that lead to the unperturbed mixing 
and argued that, in order naturally to preserve the good prediction for O12 in these schemes while 
generating an acceptably large value for ^13, it is necessary that the perturbations should leave a 
residual Z2 factor unbroken. This corresponds to the case that a single bilinear mixing term, I'a'^b 
or i^ai^c is dominant. We supported this contention by studying two examples which involved more 
general mixing and showed that, to avoid fine tuning, they reduced to single bilinear dominance. 

We also constructed a scheme that has only the atmospheric neutrino massive in the unperturbed 
case, while both the solar neutrino mass and ^13 are generated by first order mass perturbations. 
This has the merit of explaining the comparable magnitudes of the solar neutrino mass and the mass 
perturbation that generates ^13, while preserving the phenomenology of the previous schemes. 

To determine this phenomenology we developed the perturbative mixing analysis for TBM and 
for general Golden ratio schemes, and determined the resulting correlations between the magnitude 
of ^13, the other mixing angles and the Dirac CP violating phase. In doing so, guided by the mixing in 
the quark sector, we also allowed for a range of mixing in the charged lepton sector. For comparison 
we also determined the correlations for the case that ^13 comes entirely from the charged lepton 
sector with the neutrino sector being given by pure TBM, GR or BM mixing. 

For f(jZ/b or i^ai^c mixing perturbations, the correlations originate from 



to leading order in s^^. To obtain a predicted range for the CP violating observable 5 in terms of 
the measured mixing angles Oij, this relation must be corrected with terms from the charged lepton 
sector involving §23 and s\2- The accuracy of this prediction is limited by both the experimental 
uncertainties of ^23,13 and by the unknown phases that enter the contributions from the charged 
lepton sector. Using the present uncertainties on ^23,13 j the results of our analysis are summarised 
in Table [21 The differences between the predicted ranges of 6 in the various schemes largely reflect 
the different values of C in eq§^. The (TBM, GRl, GR2) schemes give C = (0.71,0.62, 1.24) for 
ab mixing and C = —(1.41, 1.62,0.81) for ac mixing. Thus the ranges for TBM and GRl are similar 
for both ab and ac mixing, while the GR2 prediction is more distinct. 

Very significant improvements are expected in the measurements of 9i2 and ^23 hi the coming 
years. In ab and ac mixing schemes, 612 is unperturbed, so that a reduced uncertainty will tell us 
whether one of these schemes, with a small contribution from 6[2, is allowed, as shown by the 
column of Table [21 and may serve to distinguish between TBM and GR schemes. On the other hand, 
a reduction in the experimental uncertainty in ^23 will lead to a significant increase in precision in 
the predicted ranges of 5, compared to the present predictions shown in Table [21 although a residual 
uncertainty from the unknown phases in the charged lepton contributions will remain. 

There are two competing ideas for understanding the large neutrino mixing angles 612 and ^23^ 
fiavor symmetries and anarchy |37l [381 139j . While the recent discovery of a relatively large value 





(52) 
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of has certainly increased the hkehhood of anarchy, the competition is far from over. Future 
precision neutrino experiments will greatly reduce the uncertainties on Oij, determine whether the 
hierarchy is normal or inverted, and finally measure CP violation, providing it is not suppressed. In 
this paper we have argued that such a program could yet uncover a very simple underlying structure 
of lepton flavor symmetry. 
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Table 2: Predictions for the allowed ranges of the Dirac CP violating phase 6 and the mixing angles. 
Two columns are shown for the prediction of 5 corresponding to the la and 3a ranges of S23 and 
S12 of the fit of Gonzales-Garcia et al [6]. For comparison, the final 2 rows show the fit to 6 and S12 
from data (c.f. Table [T|). A dash for a mixing angle indicates that the fitted value has been used to 
determine the allowed range of the other parameters. A dash for the phase indicates no solutions 
possible. 
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